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speed,  the  interplanetary  magnetic  field  (IMF)  magnitude  and  polar  angle,  we  have 
made  a  major  advance  in  predicting  both  flare-induced  storms  and  recurrent  storms. 

Furthermore,  it  is  demonstrated  that  the  prediction  scheme  ran  be  calibrated 
using  the  interplanetary  scintilation  (IPS)  observation,  when  the  solar  disturbance 
advances  about  half-way  to  the  earth. 


It  is  shovm,  however,  that  we  are  still  far  from  a  reliable  prediction  scheme. 
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Figure  Captions 

Figure  1.  Geomagnetic  and  ionospheric  storms  are  caused  by  two  types  of  solar 
phenomena.  The  first  one  is  transient  activities,  such  as  solar  flares,  prominence 
disappearances,  coronal  mass  ejections,  while  the  second  one  is  a  high  speed  stream 
from  a  coronal  hole  which  is  the  least  active  region  of  the  sun.  Thus,  the  storm 
prediction  scheme  must  take  both  phenomena  into  account. 

Figure  2.  There  is  a  considerable  spread  in  the  transit  time  (Ts)  of  solar  disturbances 
associated  with  solar  flares  to  reach  the  earth.  The  histogram  shows  that 
geomagnetic  storms  can  occur  about  20  hours  after  a  flare,  but  can  be  as  late  as  70 
hours.  Therefore,  although  the  mean  value  of  the  transit  time  is  43  hours,  such  a 
statistical  study  is  little  use  for  predicting  the  onset  time  of  a  geomagnetic  storm  after 
a  specific  flare. 

Figure  3.  The  modern  geomagnetic  storm  prediction  scheme  must  depend  on  the 
latest  progress  in  solar-terrestrial  physics,  not  on  statistical  results  shown  in  Figure  2. 
It  is  based  on  the  fact  that  a  geomagnetic  storm  results  from  increased  electric 
currents  around  the  earth,  which  are,  in  turn,  caused  by  an  increased  power  of  the 
solar  wind-maguetosphere  generator. 

Figure  4.  The  modern  geomagnetic  storm  prediction  scheme  consists  of  three  steps. 
After  parameterizing  and  quantifying  a  solar  flare,  three  vital  quantities  (V,  B,  0)  are 
computed  at  the  earth's  location  in  step  1.  In  step  2,  the  three  values  are  used  to 
compute  the  power  of  the  solar  wind-magnetosphere  generator  as  a  function  of  time 
Since  the  power  is  related  to,  empirically  or  theoretically,  to  the  geomagnetic  storm 
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indices  (AE,  Pst,  Kp),  it  is  possible  to  predict  the  the  indices  as  a  function  of  time  in 
step  3. 

Figure  5.  It  is  important  to  be  able  to  calibrate  the  storm  prediction  scheme  when 
the  responsible  solar  disturbance  reaches  about  half  way  to  the  earth.  By  monitoring 
celestrial  radio  sources  (A,  over  the  sky,  it  is  now  possible  to  infer  the 

geometry  of  the  advancing  solar  disturbance,  since  the  solar  disturbance  (the  plasma 
cloud)  causes  intense  scintillation  of  radio  sources.  Details  of  the  scheme  is  shown 
in  Figure  6. 

Figure  6.  Based  on  the  six  flare  parameters  of  a  specific  flare,  the  modern  storm 
prediction  scheme  can  determine  the  geometry  of  the  advancing  solar  disturbance. 
The  prediction  scheme  can  then  produce  the  predicted  sky  map  of  interplanetary 
scintillation  (IPS)  at  a  given  time,  say  24  hours  after  a  specific  flare.  When  the 
observed  IPS  sky  map  becomes  available,  both  sky  maps  can  be  compared.  (In  this 
particular  example,  the  IPS  was  expected  to  occur  in  the  western  sky,  but  the 
observed  high  (H)  scintillation  occurred  in  the  eastern  sky).  By  this  calibration 
method,  it  was  found  that  the  responsible  flare  was  misidentified  in  the  first 
computation;  during  an  active  period  of  sun,  many  flares  take  place  within  a  few 
hours  over  the  solar  disk.  The  correct  identification  of  the  responsible  flare  was 
confirmed  by  finding  that  another  flare  can  reproduce  the  observed  IPS  sky  map. 
Based  on  the  new  six  parameter  for  that  flare,  the  prediction  scheme  is  carried  out  to 
demonstrate  that  the  observed  variations  of  the  solar  wind  speed  is  very  close  to  the 
predicted  one. 

Figure  7/a.  Comparison  of  the  observed  solar  wind  speed  by  an  earth-bound  satellite 
and  the  predicted  one  for  the  September  25  -  October  1,  1978  storm  event.  In  spite  of 
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a  very  complicated  series  ot  events  on  the  sun,  the  three  shock  structures  are  well 
reproduced  by  the  prediction  scheme  (see  also  Figure  6). 

Figure  7B.  The  predicted  geometry  of  the  advancing  interplanetary  disturbance  for 
the  storm  event  illustrated  in  Figures  6  and  7A.  S;  Sun,  H.  earth. 

Figure  7C.  The  predicted  IPS  sky  map  of  the  advancing  interplanetary  disturbance 
for  the  storm  event  illustrated  in  Figures  6,  7A  and  7B.  Note  that  the  second  IPS, 
caused  by  a  new  flare,  began  on  September  25, 1978. 

Figure  7D.  The  observed  IPS  sky  map  of  the  advancing  interplanetary  disturbance 
for  the  storm,  event  illustrated  in  Figures  6,  7A,  7B  and  7C.  The  high  scintillation 
area  is  indicated  by  H.  This  observation  was  made  by  the  Cambridge  IPS  group. 

Figure  8A.  In  the  modern  storm  prediction  schem.e,  the  geomagnetic  storm  indices 
(AE,  Dst,  Kp)  are  estim.ated  empirically  on  the  basis  of  the  computed  power  of  the 
solar  wind-magnetosphere  generator  (Figure  4).  This  is  because  there  is  so  far  no 
reliable  theoretical  formulation  to  relate  the  power  to  the  indices.  Flowever,  here 
we  attempt  to  set  up  an  equivalent  electrical  circuit  for  the  magnetosphere.  Thus,  by 
providing  the  predicted  power  as  a  function  of  time,  it  will  be  possibb  to  determine 
the  geomagnetic  indices. 

Figure  8B.  For  a  given  voltage  variation  (OCT  (D)/  the  equivalent  circuit  can  predict 
variations  of  the  resistance  Rt  in  the  magnetotail,  the  polar  cap  potential  drop  (^ps) 
and  the  energy  (Wb)  released  in  the  polar  ionosphere.  From  these  predicted  values, 
it  is  possible  to  infer  the  geomagnetic  indices  (see  Figure  9). 
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Figure  9.  The  equivalent  circuit  in  Figure  8A  can  predict  the  potential  drop  (  ^  ps) 
across  the  polar  cap.  There  are  two  important  empirical  relationships  betwi'en  4)  ps 
and  the  AE  index,  obtained  by  satellite  observations  (Weimer,  et  al.)  and  ground- 
based  observations  (Ahn,  et  al.).  Using  this  established  relationship,  it  is  possible  to 
infer  the  AE  index  as  a  function  of  time  from  (t>  ps- 

Figure  10.  We  can  advance  the  present  geomagnetic  storm  prediction  scheme  to 
include  the  prediction  of  ionospheric  storms.  Once  the  power  of  the  solar  wind- 
magnetosphere  generator  is  predicted  in  step  2  (Figure  4),  it  is  possible  to  predict  the 
potential  drop  ^  pc  across  the  polar  cap,  as  illustrated  in  Figures  8A,  8B  and  9.  Our 
prediction  scheme  can  also  predict  the  size  of  the  auroral  oval.  Thus,  for  a  given 
date  and  time,  the  ionosphere  storm  code  can  predict  the  flow  of  ionospheric  plasma 
across  the  polar  region  and  the  electron  density  profile  at  a  desired  time.  Radio  ray 
paths  can  be  computed  on  the  basis  of  the  predicted  electron  density  profile.  In  this 
case,  a  high  ‘h  pc  caused  an  intense  flow  of  ionospheric  plasma  from  the  sunlit 
hemisphere  to  the  dark  hemisphere  across  the  polar  cap. 

Figure  11.  Example  of  the  computed  electron  density  profile  at  an  altitude  of  300 
km.  There  is  a  high  density  in  the  dayside  ionosphere  (caused  by  the  solar  radiation) 
and  also  along  the  auroral  oval  (caused  by  auroral  particle  precipitation).  The 
dayside  plasma  flows  into  the  polar  cap  by  the  dawn-dusk  potential  drop  ( ^  pc)  can 
be  seen. 

Figure  12A.  Simulation  result  demonstrating  how  a  particular  type  of  disturbance 
caused  at  a  confined  location  in  a  magnetic  flux  tube  propagates.  Such  a  disturbance 
may  be  identified  as  a  "magnetic  cloud.” 
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Figure  12B.  Our  prediction  scheme  has  bet-p.  improved  to  ihe  point  to  be  ahl.'  to 
predict  the  solar  wind  speed  V  and  ihe  interplanetar v  p'.: .’.p':;,-  PelP  (IMF.) 
magnitude  3  for  simple  events.  The  simplest  v-vent  is  illusinti. d  i.i  a;  in.  v.  hich  'he 
advancing  shock  wa^'ets)  compresses  the  existn-.'i  inier[>lari  jar'-  m.agrulic  field  in 
tne  vicinity  of  tlie  earth.  !f  the  IMF  field  is  di'-ected  south vvard  and  iias  a  magnitude 
of  5  nT,  'he  arrival  of  ihe  shock  wave  can  enmpres,’.  tee  fiei'l,  eev'sing  e  large  (~  15 
nT)  souihward  pointing  field.  Such  a  high  southwani  iie'i.i  car  rertaivilv  cause  an 
intense  geomagnetic  storm.  Ftowever,  the  advancing  disturbances  are  often  .much 
more  complicated.  In  an  exactly  opoposite  case,  an  enhanced  northward  directed 
field  will  .r'i'cluce  geon'iagnetic  activity  A  number  of  possibilities  ..ire  shown  here. 

Figure  13.  During  1974,  there  were  tv/o  large  coronal  halos,  from  which  two  high 
speed  streams  em,anated.  During  a  single  solar  rotation  {27  days),  two  peaks  of  the 
solar  wind  speed  were  seen  (lop).  The  next  two  rows  sni.w';-,  th.e  IMF  magnitude 
and  the  two  angles  of  the  IMF  tthe  latitude  angle  theta  and  the  azimutnal  angle  phi). 
On  the  basis  of  these  observed  quantity,  the  power  {  r  )  c^f  the  solar  wind- 
magnetosphicre  generator  is  computed.  The  lust  two  r.vv-.,  .how  the  AF  and  Dst 
index,  respectivity.  There  is  a  f.ainy  good  similarity  betvvi  cn  c  and  AE,  suggesting 
that  t  can  predict  tlie  envelope  of  the  AE  index. 

Figure  14A.  The  prediction  of  a  high  speed  stream  depends  cnticaiiy  on  the 
magnetic  field  distribution  on  the  so-called  "source  surtaee,"  a  spherical  surface  of 
radius  of  2,5  solar  radii.  In  the  right  colum.n,  sunspot  number  betw'een  1976  and 
'986  IS  showm.  In  the  middle  column,  .shows  how  the  iOagnotiC  neutral  line  (or  the 
magnetic  equator'  on  !he  source  .siuface  changed  during  dm  sunspot  cycle.  The  left 
column  shows  the  polarity  of  tl\e  magnetic  field  and  the  magnetic  equator  in  solar 
longitude-latitude  map  (the  magnetic  field  is  directed  im.s'ard  in  ihe  hatched  area). 
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Figure  14B.  In  order  to  predict  the  geometry  of  the  magnetic  neutral  line,  a  simple 
method  is  developed  to  represent  the  magnetic  neutral  line.  The  method  uses  a 
central,  axially  aligned  dipole  and  two  (or  at  most  four)  dipoles  on  the  photosphere. 
The  observed  and  reproduced  magnetic  neutral  lines  are  compared  for  C  arrington 
rotation  1666. 

Figure  15.  Once  the  magnetic  neutral  line  is  reasonably  reproduced,  the  prediction 
of  the  magnetic  neutral  line  is  reduced  to  the  prediction  of  changes  of  the  two 
photospheric  dipoles  in  Figure  14B.  The  figure  shows  the  observed  and  reproduced 
magnetic  neutral  lines  for  Carrington  rotations  1661-1664.  In  these  cases,  four 
photospheric  dipoles  are  used. 

Figure  16.  In  this  figure,  we  follow  the  location  (latitude,  longitude),  the  azimuth 
angle  and  magnitude  of  the  photospheric  dipole  no.  3  in  Figure  15.  One  can  see  that 
the  changes  of  the  dipole  were  not  erratic  Thus,  it  was  possible  to  extrapolate 
changes  of  the  dipole  quantities  and  thus  predict  changes  of  the  magnetic  neutral 
time  during  the  successive  Carrington  rotations. 
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Development  of  h  Numerical  Scheme  m  Predict  Geomagnetic  Stia  ms 
After  intense  Solar  Events  and  (Ejomr.gnetic  Arti’/itv  'Z7  Days  in  Advan;  t 

S,-I.  Akasotu 

i epiivs-r.'.-i  irisi';. 

L'liiver-'iiy  of  Aiaska 
:\-.irhouk'  Ad.  P977r 

SUMM/.AV 

Tiie  (  onlrac!  nas  two  major  ohjectives.  The  lirct  (ibn-etive  is  to  (>sia*Mish  th' 
computer  schojao  s-, nich  can  pred’ci  the  occurrence  and  deveiopnieat  of  a  lajor 
geomagnetic  storm  after  an  '  .tense  solar  i.are.  Since  the  intensity  of 
geo:-';;^  nc..;c  ■forn:  is  given  in  terms  of  the  gcmmagnetic  indices  such  as  AE,  Kp 
ana  D.^n  a  ....  .ve  storm  prediction  mast  be  able  to  provide  the  geomagnetic 

indices  as  a  tunciion  of  time  for  a  given  solar  flare. 

The  pr  edic  tion  .sch.C'me  consists  of  three  steps.  The  first  step  is  an  attempt 
to  compute  th.  srd  .'-  wuid  speed  :V)  the  inter planct-ry  um-r  i  -lie  field  magnitude 
and  llu'  polar  angle  (b;.  'I’he  second  step  computes  the  pmver  (Pi  generated  by 
the  solar  wnd noagrie*  espio  rc  hiteraition.  !n  the  ‘hir'd  steo  tsio  empiricai 
relationsnip  between  P  aitd  the  gcomagnetiv:  indices  i-.  used  to  eonipuio  the 
geomagnetic  indices  as  a  function  of  time.  At  the  present  tur.e,  it  is  not  pc'^sible  to 
predicl  the  polar  ajigle  cf  6  of  the  interrilanelary  maeuc.ie  hoih.  This  reduces  our 
scheme  to  predict  t’nc  maximum  possible  power  Fni. 

•  We  have  .ompieled  ihc'  solar  wind  cod.e  which  p.-cdicto  the  sober  wind  speed 
V  and  the  interplanetary  magnetic  Held  magnitude  B  at  the  front  of  the 
magnetospliere  even  for  a  complicated  situation  in  which  many  flares  occur  in  a 
relatively  short  period.  The  two  quantities  aliovv  us  to  compute  the  maximum 
posshdf'  pic,'  t  r  B,  ;  (h  th(  .njiar  c^nnd  mag.mtospiier'  gencronoi' 

•  In  the  prediction  scheme,  it  is  important  to  he  able  to  compare  the  observed 
and  computed  geo  iiet.y  of  the  shock  wa  rn  at  a  midpoint  hetwi'cn  fhe  sun  and  the 


earth  We  have  found  that  the  interedaiietary  scint iliat -ju  Mr;-  ^  m  fwr 

such  a  purpose.  Thus,  we  hav»  successfully  incorporated  results  from  tlic 
int'^'rplanetary  scintillation  (1PS»  in  the  prediction  scheme.  In  the  cod(‘.  ve  can 
calibrate  the  initial  conditions  of  a  soiaz’  flare  using  the  IPS  observat  ions 

•  Since  we  can  compute  P^tt)  as  a  function  of  time  and  .-in.-e  ilior^-  no 
empirical  relations  between  Pit)  and  AE, Kp/Dst.  is  possible  t(j  predict  the 
maximum  AE/Kp/li)st  a&  a  fnuction  of  time. 

'  Since  there  is  ar  err.pirical  relationship  between  P  and  the  cross  polar  cap 
potential  ^p,-,  it  is  also  possible  to  infer  time  variation.'  of  the  ntaximum  s  ;duo  of 
'h/)c  for  a  given  time  vanation  of  Pn^ 

.  We  have  also  developed  an  equivalent  circint  for  the  magnetosphere.  For  a 
given  input  function  Pit)  or  Pni(t),  it  is  possible-  to  predict  AE.fKp.fDst  (or  the 
maximum  of  AE/KP/Dst)  as  a  function  of  time. 

•  We  have  also  developed  a.  code  to  predict  the  distribution  of  tne  F-laye.*' 
electron  density.  The  F-layer  electron  density  in  the  polar  cap  depends  crucially 
on  the  potential  dnzp  across  the  polar  cap. 

The  second  objective  of  *his  project  is  to  advance  the  prediction  of  27-day 
recurrent  storms. 

•  We  have  found  that  this  project  is  reduced  to  predict  both  the  solar  wind 
speed  and  the  I?4F  magnitude  27  days  in  advance  and  that  this  task  i.s  further 
reduced  to  predict  the  geometry  of  the  .'■eulral  line  cm  ..he  solar  sc  lire. surf.ace  hi 
sphe»*ioal  surface  of  2  5  solar  radii). 

•  We  have  devised  a  multi-dipole  method  to  reproduce  the  neutral  line 

•  Thus,  the  prediclucu  is  finally  reduced  lo  prod’ct  the  growth  and  dr-cay  of 
the:  dipoles  in  low  latitude  of  the  phetosphere  We  have  made-  iouch  uregre  ss  in 
under.standing  the  nature  of  these  dipide.s. 


In  summary,  altogether  we  have  developed  a  computer  code  to  execute 
what  is  described  in  the  above,  our  present  prediction  scheme  may  bo  considered 
to  be  a  skeleton  structure'  for  the  future  scheme.  We  believe  that  most  of  the 
elements  in  the  prediction  scheme  are  assembled  and  arc  properly  connected. 
The  method  requires  much  improvement  in  the  future,  which  depends  on  the 
progi'css,  both  observational  and  theoretical,  in  the  field  of  solar-terrestrial 
physics.  These  are  summarized  m  the  last  section  "Concluding  Remarks." 


1.  INTRODUCTION 
1.1  Genera!  Review 


Solar-terrestrial  physics  has  been  advanced  considerably  during  the  last 
two  decades.  The  progress  has  been  prominent  in  almost  every  discipline, 
namely  studies  of  solar  flares,  interplanetary  disturbances  and  magnetospheric 
disturbances.  A  transient  phenomenon  called  "coronal  mass  ejection  (CME)”  has 
been  discovered  (cf.  Gosling,  et  al.,  1974;  Hildner,  et  al.,  1975,  1976;  Sheeley,  et  al., 
1980;  Wagner,  1984)  and  its  association  with  soft  x-ray  events  and  other 
phenomena  have  been  investigated  (Sheeley  et  ai.,  1983,  1984;  Kahler  et  al.,  1984a, 
b;  Burlaga  et  al.,  1982).  The  source  region  of  high  speed  solar  wind  streams, 
namely  the  "coronal  holes,"  have  also  been  identified  (cf.  Zirker,  1977).  A  study  of 
interplanetary  disturbances  has  been  greatly  advanced  by  the  availability  of  in  situ 
observations  of  the  solar  wind  and  the  interplanetary  magnetic  field  (IMF)  by  deep 
space  probes  and  earth-bound  satellites.  The  interaction  between  the  quiet  time’ 
solar  wind  and  a  high  speed  solar  wind  stream,  namely  the  stream-stream 
interaction,  has  been  studied  extensively,  both  observationally  and  theoretically 
(cf.  Dryer  and  Steinolfson,  1976;  D'Uston,  et  al.,  1981). 

In  addition  to  such  progress  in  each  discipline  of  solar-terrestrial  physics, 
there  have  been  considerable  efforts  in  finding  one-to-one  relationships  between 
solar  activity  and  interplanetary  disturbances  and  also  between  interplanetary 
disturbances  and  magnetospheric  substorms  and  storms.  An  extensive  study  of 
the  relationship  between  the  large-scale  magnetic  field  in  the  photospheric  level 
and  the  interplanetary  magnetic  field  in  the  photospheric  level  and  the 
interplanetary  magnetic  field  has  been  made  (cf.  Wilcox  et  al.,  1980;  Hoeksema  et 
al.,  1982,  1983),  The  relationship  between  solar  activity  (such  as  flares  and  CMEs) 
and  shock  waves  in  interplanetary  space  has  been  investigated  (Chao  and 
Lepping,  1974;  Joselyn  and  Bryson,  1980;  Sheeley  et  al.,  1985).  The  so-called 
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"Sudden  Disappearing  Filaments  (SDFs)"  have  been  added  as  a  new  possible 
source  of  interplanetary  magnetospheric  disturbances  (Schwenn  et  al.,  1983; 
Joseljm  and  McIntosh,  1981;  Wright  and  McNamara,  1983;  Sanahuja  ct  al.,  1983). 
In  solar  wind-magnetosphere  interaction  studies,  some  of  the  key  physical 
quantities  in  the  energy  transfer  processes  have  been  identified.  Among  these, 
the  north-south  component  of  the  IMF  is  found  to  play  a  crucial  role  in  the  energy 
transfer  (cf.  Arnoldy,  1981;  Tsurutani  and  Meng,  1972;  Russell  and  McPherron, 
1973;  Meng  et  al.,  1973;  Clauer  et  al.,  1981;  Akasofu,  1981;  Baker  et  al.,  1981; 
Meloni  et  al.,  1982). 

In  spite  of  such  progresses  in  each  discipline  and  in  some  interdisciplinary 
areas,  our  understanding  of  some  of  the  most  crucial  aspects  of  solar-terrestrial 
physics  has  not  been  very  much  improved.  One  of  such  examples  is  the  so-called 
"driver-gas"  which  is  supposed  to  be  ejected  during  a  solar  event  and  to  cause  an 
interplanetary  shock  wave  which  causes,  in  turn,  a  geomagnetic  storm. 
However,  the  identification  of  the  "driver-gas"  in  interplanetary  space  has  not 
been  very  conclusive.  First  of  all,  only  a  few  plasma  ejectas  identified  in  the  Ha 
photographs  have  been  observed  to  leave  the  field  of  view  of  the  camera;  some  of 
them  simply  rain  back  down  to  the  chromosphere.  Tracing  of  type  IV  radio 
bursts  which  are  supposed  tc  be  generated  by  energetic  electrons  in  the  'ejecta' 
have  been  limited  to  about  several  solar  radii.  Except  for  a  fairly  high  correlation, 
the  physical  relationship  between  a  CME  and  a  shock  wave  has  not  been 
established  yet.  Is  a  CME  an  early  stage  of  the  shock  wave,  a  driven  gas  or 
neither  of  them? 

There  has  also  been  no  agreed  common  signature  of  the  driver  gas  from 
spacecraft  measurements.  Some  workers  consider  that  a  high  ratio  of  He/H  (or  a 
high  concentration  of  helium)  after  the  passage  of  the  shock  waves  is  an 
important  signature  of  the  "driver  gas"  (Hirshberg  et  al.,  1970,  1972a,  b;  Gosling  et 
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al.,  1980;  Borrini  et  al.,  1982a,  b;  Zwickl  et  al.,  1983),  Others  consider  a  prominent 
decrease  of  the  proton  temperature  (Gosling  et  al.,  1973;  Zwickl  et  al.,  1983),  or  a 
very  steady,  large  magnitude  of  the  interplanetary  magnetic  field  (Smith,  1983)  as 
the  most  important  signature.  The  so-called  "magnetic  cloud”  proposed  by  Klein 
and  Burlaga  (1982)  may  also  be  added  as  another  proposed  signature.  Bi¬ 
directional  streaming  of  solar  wind  electrons  has  also  been  considered  as  an 
evidence  for  a  closed  field  structure  within  the  drive  gas  (Bame  et  al.,  1981).  For  a 
recent  review  on  the  plasma  properties  of  the  solar  wind,  see  Neugebauer  (1983). 

A  number  of  MHD  simulation  studies  have  also  been  conducted  to  predict 
how  the  density,  temperature  and  velocity  of  the  disturbed  solar  wind  and  IMF 
vary  behind  the  interplanetary  shock  waves.  The  common  practice  in  this 
particular  simulation  study  is  to  allow  the  duration  of  heat  input  during  a  flare  to 
be  a  variable.  An  impulsive  heat  input  of  duration  of  30  min  -  1.5  hours  into  the 
solar  atmosphere  simulates  the  generation  of  a  blast  wave,  while  a  constant  heat 
input  lasting  for  at  least  40  hours  simulates  a  continuous  'piston-like'  process 
with  a  shock  front  (cf.  Hundhausen  and  Gentry,  1969a,  b;  Dryer,  1970).  The 
simulation  studies  of  the  blast  waves  indicate  that  it  is  followed  by  a  large 
depression  of  the  density  and  of  the  IMF  magnitude  (cf  Wu  et  al.,  1976;  D’Uston  et 
al.,  1981;  Wu  et  al.,  1983;  Dryer  et  al.,  1984). 

However,  such  a  prominent  depression  of  both  the  proton  density  and  the 
IMF  magnitude  after  the  passage  of  the  shock  waves  is  not  a  common  occurrence, 
although  the  proton  density  can  drop  after  the  passage  of  several  major  shocks 
(Sheeley  et  al.,  1985).  For  this  and  other  reasons,  Borrini  et  al.,  (1982b)  concluded 
that  there  is  little  evidence  that  the  blast  waves  actually  exist  in  interplanetary 
.space  and  that  the  occurrence  of  mass  ejection  events  (Hildner  et  al.,  1976)  is  more 
than  sufficient  to  account  for  all  shocks  observed  at  1  AU.  Recently,  Smart  and 
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Shea  (1985)  suggested  that  most  of  the  shock  waves  are  initially  driven  to  a 
distance  of  0.1  AU  and  then  propagate  as  a  blast  wave. 

Recently,  Zwickl  et  al.  (1983)  showed  that  of  54  shocks  observed  from  August 
1978  to  February  1980,  9  events  were  followed  by  a  clearly  identifiable  decrease  of 
the  temperature  and  an  increase  of  the  He/H  ratio,  although  the  source  activity 
for  those  chosen  9  events  were  not  identified.  Zwickl  et  al.  (1983)  also  added  an 
increase  of  the  IMF  magnitude  as  a  signature  of  the  driver  gas.  Borrini  et  al. 
(1982b)  inferred  that  a  high  ratio  of  He/H  arises  from  flares  in  the  vicinity  of  the 
central  meridian,  although  their  paper  does  not  show  any  supporting  evidence  for 
such  an  inference. 

Geomagnetic/ionospheric  storms  can  be  classified  into  two  types, 
depending  on  their  origin  in  the  solar  atmosphere  (Fig.  1).  The  first  type  is 
associated  with  intense  eruptions  in  the  solar  atmosphere,  which  are  generally 
called  solar  flares.  It  is  known,  however,  that  not  all  flares  cause 
geomagnetic/ionospheric  storms.  Recent  studies  have  revealed  that  solar  flares 
followed  by  soft  x-ray  emissions  of  a  long  duration  (of  several  hours)  cause 
geomagnetic/ionospheric  storms  (Sheeley  et  al.,  1983,  1984).  It  is  expected  that 
such  a  flare  ejects  a  gas  cloud  into  interplanetary  space.  As  the  driver  gas 
advances  into  interplanetary  space,  it  generates  an  interplanetary  shock  wave.  A 
geomagnetic/ionospheric  storm  begins  when  the  shock  wave  and  the  driver  gas 
collide  with  the  earth’s  magnetic  field  The  development  of  a  geomagnetic  storm 
depends  on  several  physical  parameters  of  the  driver  gas  and  its  relative  path 
with  respect  to  the  earth.  As  we  reviewed  earlier,  our  knowledge  on  the  driver 
gas  is  still  too  poor  to  be  used  in  predicting  the  development  of  a  geomagnetic 
storm. 

The  second  type  of  geomagnetic/ionospheric  storms  is  associated  with  an 
intense  solar  flare  wind  flow  from  very  quiet  regions  of  the  solar  corona.  Such 
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regions  look  very  dark  when  the  corona  is  observed  by  a  soft  x-ra>  camera.  Ii  is 
for  this  reason  that  these  regions  are  called  'coronal  holes.'  It  is  not  known  at  all 
why  the  most  quiet  region  of  the  solar  corona  can  produce  the  most  intense  solar 
wind.  Hakameda  (1987)  and  Wang  and  Sheeley  (1990)  showed  that  the  highest 
speed  wind  may  be  related  to  very  limited  regions  on  the  photosphere.  The  life 
time  of  coronal  holes  varies  considerably,  from  less  than  one  month  to  more  than 
12  months.  A  long  lasting  coronal  hole  causes  geomagnetic  storms  with  an 
interval  of  approximately  27  days,  as  the  sun  rotateo  once  in  25  days  and  the  earth 
revolves  around  the  sun  in  the  same  direction. 

1.2  Needs  for  a  New  Prediction  Scheme 

In  the  past,  most  prediction  effort  of  geomagnetic/ionospheric  storms  had 
been  based  on  statistical  results.  For  example,  the  prediction  of  the  onset  time  of  a 
geomagnetic  storm  is  made  on  the  basis  of  a  statistical  result  that  geomagnetic 
storms  begin  most  frequently  about  43  hours  after  a  responsible  flare.  Figure  2 
shows  the  histogram  which  shows  the  time  interval  between  flare  onset  and 
storm  onset.  Although  such  a  histogram  is  useful,  the  spread  is  too  large  to  be 
useful  in  predicting  the  onset  time  of  a  geomagnetic  storm  after  a  specific  solar 
event. 

Like  the  modern  weather  forecasting,  the  modern  prediction  scheme 
should  be  based  on  numerical  methods  for  individual  events  and  is  based  on  the 
following  facts  (Fig.  3): 

(1)  The  interaction  between  the  solar  wind  and  the  earth's  magnetic  field 
constitutes  a  generator. 

(2)  The  power  of  the  generator  depends  on  the  solar  wind  V,  the  magnitude  B 
and  the  polar  angle  9  of  the  interplanetary  magnetic  field  (Perreault  and 
Akasofu,  1978;  Akasofu,  1981). 
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GEOMAGNETIC  AND  IONOSPHERIC  STORMS 
First  Type  -  associated  with  Solar  Flares 
Second  Type  --  associated  with  Coronal  Holes 


b 


FIGURE  1 
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MODERN  GEOMAGNETIC-IONOSPHERIC  STORM 

PREDICTION  SCHEME 
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igure 


P  (Mwatts)  =  20  V  (km/sec)  B^(nT)  sin"^  (6/2) 

(3)  Thus,  having  a  higher  speed  V  and  greater  magnitude  B  of  the 
interplanetary  magnetic  field,  a  'gusty'  solar  wind  (the  shock  wave/driver 
gas  or  of  the  solar  wind  flow  from  the  coronal  hole)  tends  to  increase  the 
power  of  the  solar  wing-magnetosphere  generator. 

(4)  The  resulting  intensified  electric  currents  generate  complex  magnetic  fields 
which  are  the  cause  of  a  geomagnetic  storm. 

(5)  Part  of  electrons  carr3dng  the  intensified  currents  collide  with  upper 
atmospheric  atoms  and  molecules  and  are  responsible  for  causing 
ionizations  there,  causing  auroral  phenomena. 

(6)  An  increased  potential  drop  associated  with  an  enhanced  power  drives  the  F- 
layer  ionization  in  the  dayside  to  drift  into  the  polar  cap,  causing  major 
changes  in  the  distribution  of  the  ionization. 

A  geomagnetic/ionospheric  storm  and  an  auroral  display  are  nothing  but 
different  manifestations  of  the  intensified  currents.  The  higher  the  power  of  the 
generator  is,  the  higher  is  the  current  intensity,  causing  more  intense 
geomagnetic  disturbances  and  brighter  auroral  displays.  This  is  why  both  a  very 
intense  geomagnetic/ionospheric  storm  and  a  great  auroral  display  occur 
together. 

Therefore,  the  prediction  of  geomagnetic/ionospheric  storms  and  the 
aurora  is  reduced  to  predicting  the  power  of  the  generator  which  depends  on  the 
solar  wind  speed  V,  the  solar  wind  magnetic  field  B  and  its  polar  angle  0.  The 
first  task  of  this  project  is  thus  to  develop  a  method  to  predict  the  above  three 
quantities  (V.  B.  9)  for  the  interplanetary  shock  wave/driver  gas  and  for  an 
intense  solar  wind  flow  from  a  coronal  hole. 
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In  section  2,  we  de  cnho  the  p’'ojrress  we  have  made  in  predicting  major 
geomagnetic  storms,  in  particular  our  efforts  to  cooperate  the  interplanetary 
scintillation  HPS'  ouh.  ■  va a  .n.  As  wt  shaH  describe  in  more  detail  later,  the  IPS 
observation  servos  t’t  o  pu;  of  checking  the  prediction  at  a  midpoint  between 
the  sun  and  iho  earth  i.,  seirir;;  3,  we  describe  the  progress  w'e  have  made  in 
predicting  the  2'' -day  r  current  geomagnetic  storms  We  show  that  this 
particular  project  is  lediued  to  predict  the  geometry  of  the  neutral  line  on  the 
source  surface  of  ihe  sua  on  the  basis  of  the  observed  neutral  lines  during  the 
past  Carrington  totatnin.  Thus,  much  of  our  effort  has  been  concentrated  in 
understanding  the  iiu  gnetic  holds  on  the  source  surface  and  their  variations. 


The  modern  prodicnon  schenoe  should  be  able  to  provide  not  only  just  the 
occurrence  and  the  maximum  intensity  for  a  given  storm,  but  also  the  time 
development  of  th;.  .'tonn  Since  the  intensii.y  of  a  geomagnetic  storm  field  at  a 
given  time  can  be  expressed  in  terms  of  the  geomagnetic  indices,  we  must  be  able 
to  predict  the  geomt  gnotif  indices  as  a  function  of  time  (such  as  AE,  Dst,  Kp, 
etej.  namely  AC  '.).  D-ihC  K.ii ;  etc.  Toward  this  goal,  we  have  developed  a 


three-step  process;  'see  Akssofu  and  Fry  (1986)  and  Figure  4). 

2.1  STEP  1 

STEP  1:  On  the  basis  if  the  ob.scrvation  of  a  particular  solar  flare,  we 

determine  ttic  iolh.wing  parameters; 

1.  d  ime  '  UT),  dal  e  of  the  occurrence 


-V  ),  lairiudo  (tNkS)  of  flare 


3.  Initial  speed  of  the  solar  wdnd’' 

4.  Area  covered  by  the  flare 


5.  Duration 


The  quantities  with  *  depend  on  the  intensity  of  a  flare  and  are  f^ivcn  in 
Table  1. 

Table  1:  Conversion  From  Flare  Importance  and 
Brilliance  to  V/r  and  Gp 


Importance/Brilliance 

(km  s  U 

Op  (*-) 

OF 

200 

20 

ON 

300 

30 

OB 

400 

40 

IF 

400 

30 

IN 

600 

40 

IB 

800 

50 

2F 

500 

30 

2N 

800 

50 

2B 

1200 

70 

3F 

700 

30 

3N 

1100 

50 

3B 

ibOO 

70 

4F 

800 

40 

4N 

1200 

60 

4B 

1600 

80 

Using  the  above  solar  flare  parameters  as  the  input,  a  computer  program 


has  been  developed  to  predict  the  solar  wind  speed  (V)  and  the  magnitude  (B)  of 
the  interplanetary  magnetic  field  at  the  front  of  the  magnetosphere  as  a  function 
of  time,  namely  V(t),  and  B(t).  When  a  number  of  flares  take  place  during  a  given 
period,  they  can  be  inputed  successively  into  the  program.  This  is  often  the  case, 
when  there  are  one  or  more  very  active  regions  on  the  solar  disk. 

STEP  2;  The  predicted  values  in  the  above  enable  us  to  estimate  the  total 
maximum  power  generated  by  the  magnetosphcric  interaction  aii 

a  . function  of  time, 
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’’’m  (Mv.'atts;  =  20  V  (km/sec)  B^(nT) 

Note  that  Pn\  ciit'feis  from  P  hy  eliminatir^  Si«'*(8/2j£' ’  0 
STEPS:  On  the  be  is  vi  'le  eeipirical  relationships  between  P  and  AE.  P  and 

Dst  P  anc  otr  ,  v,e  predict  the  geomagnetic  indices  AE(t),  Dst(t) 
and  Kp!t)  a."  a  lunition  of  time.  There  is  also  an  empirical 
relatioTisriip  '  ,  :  sve:n  ?  and  the  cro  -r-polar  cap  potential  d>pc- 

•  Using  AiPi),  (sat)  and  Kp(t)  thus  determined,  some  other 
importaia  paraaii. ters  can  a;si.  be  predicted; 

•  i  he  aianw,  Ler  (d  the  auroral  oval  can  also  be  estimated  from  the 
[iredi.d  'i  Pb(*)  and  Dst(t),  since  an  empirical  relationshi^- 
is  t  Men  ih.  f'va!  diam<3l':  end  AE/lJst  has  been  established. 

•  The  distnb’  tion  of  the  electron  density  in  the  polar  ionosphere 
alter  slorn'i  onset  can  also  be  predicted  as  a  function  of  time  by 
ktK'Witij.,  fhp  dia'Tivtcr  of  the  auroral  oval  and  the  cross-polar  cap 
pot  on  dal  C'  . 

We  have  tc:  ted  the  proftram  extensively  by  'postdicting'  many  past  major 
solar  terroGtr'il  '>v.,nts  d'inng  the  last  ten  years  or  so  and  calibrated  the 
computer  code  icf.  Akasofu  :n.d  '.^  e,  1988,  1989,  1990). 

One  of  the  diirifulde;  vv*  e-,  countered  in  applying  our  prediction  scheme  is 
that  a  niajr  r  acivu'  ivgaiUi  on  die  sun  produces  many  flares,  more  than  10  in  a 
two-week  penod  euare  a-  i  all  Dares  produce  interplanetary  disturbances,  it  has 
been  difficult  Ui  rvi  the  far  cation  alone  to  determine  which  flares  are  likely 

to  ha\e  caused  inieip!;,-:-  :  -ry  dicturh.ances. 
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FIGURE  4 


It  is  very  fortunate  that  an  analysis  of  the  interplanetary  scintillation  (IPS) 
has  now  become  availahlo.  In  the  IPS  study,  radio  waves  from  a  predetermined 
set  of  radio  galaxies  .ea  <010:  ;. rs  are  inonitoied  by  a  fixed  antenna  system  or 
systems.  Thus,  i!  a  p’osaia  i  cud  ejected  by  a  solar  flare  is  located  between  a 
radio  galaxy  and  the  n’  t!,,  a  -."ojig  scintillation  of  the  radio  wave  from  the  radio 
galaxy  is  observed  iFigorv  I',  The  fixed  antenna  system  can  scan  the  entire  sky, 
as  the  sky  rotates'  around  the  earth.  Since  the  plasma  cloud  covers  only  a 
portion  of  the  sky,  the  IPS  region  occupies  only  a  limited  region  in  the  sky  (the  sky 
mapj.  Therefore  we  have  modified  our  computer  code  to  produce  predicted  IPS 
sky  maps  which  indicate  portions  of  the  sky  where  IPS  is  expected  to  occur  for  a 
given  Hare. 


In  thus  way.  we  can  compaie  the  observed  sky  map  and  the  predicted  sky 
map  during  the  tran.  it  of  the  nlasma  cloud  from  the  sun  to  the  earth.  If  there  is 
any  major  disayrc*  ment  bc*%vcen  them,  wo  can  modify  the  input  Hare 
parameters  uni.i'  b'hh  .agree.  The  steps  to  be  taken  in  this  scheme  are 
graphically  illustr .ilc  i  i.i  Figure  fS. 

'1;  Deteimine  liiO  .-or.  o  1:  air.'Ters  from  the  observation  of  a  flare. 

(lU  Using  the  snl.i r  v;ii;d  code,  determine  the  geometry  of  the  shock  front  at 

piedeteri.ui  ‘vi  1  hi;  -  afto.  ihc  flare.  The  predeter  nined  hour  depends  on 
the  hvata  a':;)  cl  Ihc  IPS  observation. 


(iii)  Project  the  c  nn.iutec  .feiimciry  fd’the  shock  Vvave  in  che  sky  map. 
five  (lomp.jie  the  sky  ni;.  nus  produced  v/ith  the  observed  IPS  sky  map. 
tv)  If  both  <ic  v  r  -ai;.:  ogrea,  extend  the  comput'Tion  on  the  basis  of  the  same  six- 
parameter, s  are'  compote  V  and  H  and  then  power  P. 

(vii  If  the  two  :  ky  riau-  are  sm •.ita.'itialiy  diferent,  change  the  six  parameters 
or  ch(,')',e  a  (iifiVrenl  fiare. 

(vii)  Repeat  the  process  tea)  until  the  two  maps  agree  ri'asonably  well. 


(viii)  On  the  basis  of  the  six  parameters  thus  determined  finally,  compute  V  and 

B  and  then  Pm- 

Thus,  it  is  our  finding  that  the  IPS  observation  can  provide  us  with  a 
valuable  mid-point  calibration  in  the  storm  prediction  scheme.  The  IPS 
observation  is  found  to  be  far  more  useful  than  a  single  space  probe  because  we 
can  infer  a  3-D  geometry  of  the  expanding  solar  plasma  cloud. 

Since  the  present  IPS  observation  uses  a  fixed  antenna  system,  the  sky  can 
be  'scanned'  only  once  a  day  at  a  station  (using  the  earth's  rotation).  At  the 
present  time,  IPS  observations  are  conducted  at  two  locations,  Cambridge 
(England)  and  Nagoya  (Japan).  Thus,  we  can  make  the  IPS  calibration'  about 
twice  during  the  transit  of  the  solar  plasma  cloud  from  the  sun  to  the  earth.  Our 
prediction  scheme  has  been  developed  to  be  able  to  produce  expected  IPS  sky  maps 
at  any  given  time  at  both  locations.  It  is  learned  that  India  will  also  soon  join  in 
the  IPS  observation. 

We  have  tested  the  geomagnetic  storm  prediction  scheme  by  including  the 
IPS  observations  and  found  that  the  IPS  observation  can  increase  considerably  the 
accuracy  of  the  prediction  (Akasofu  and  Lee,  1989,  1990). 

In  order  to  test  the  usefulness  of  IPS  observations,  we  examined  possible 
shock  waves  (in  terms  of  time  variations  of  the  solar  wind  speed,  the  arrival  time, 
etc.)  which  should  have  been  caused  by  some  of  the  seven  flares  during  the  period 
between  September  21  and  27,  1978,  by  taking  into  account  their  location  on  the 
solar  disk  and  their  intensity.  The  coronal  hole  was  also  modeled  carefully  in 
order  to  account  for  the  observed  variations  of  the  solar  wind  speed  for  the  period 
between  22  Septernbe?  and  6  October.  Since  these  interplanetary  events  tend  to 
superpofse  upon  each  other,  the  speed,  arrival  time  and  intensity  of  a  particular 
shock  event  cannot  be  investigated  without  considering  all  responsible  flares  and 
high  speed  streams.  A  quantitative  method  is  needed  to  identify  solar  flares 
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which  are  responsible  for  the  observed  interplanetary  events.  The  method  we 
have  developed  is  capable  of  handling  several  solar  events. 

It  is  aftei  a  large  number  of  trials  and  errors  that  we  have  tentatively 
identified  the  cause-eiTecl  relationship  between  the  flares  and  the  interplanetary 
events  during  the  period  of  our  concern  (Table  2).  We  have  identified  FI  on  21 
September  as  the  cause  cf  the  shock  wave  which  was  observed  on  25  September  at 
the  Earth.  However,  this  identification  is  not  in  agreement  with  that  made  by 
Cane  et  al.  (1982).  Thus,  based  on  their  conclusion,  we  repeated  the  modeling  by 
assuming  that  F2  on  23  September  is  the  cause  of  the  shock  wave  on  25  September. 
The  former  modeling  is  called  here  Case  1  and  the  latter  Case  2  (see  Table  2). 
Table  3  shows  the  six  parameters  for  the  chosen  flares,  FI  (Case  1),  F2  (Case  2),  F3 
and  F6.  It  is  assumed  in  this  study  that  the  background  solar  wind  stream 
without  the  high  speed  ^t^eam  is  350  km/sec.  This  is  a  reasonable  choice  based  on 
the  observations  during  this  period. 
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Table  2:  List  of  Solar  Flares  Which  Occurred  Between 
21  and  27  September  1978:  Among  Them  FI,  F3  and  F6  are 
Identified  to  be  Responsible  for  SI,  S2  and  S3,  Respectively. 


Table  3:  Parameters  Employed  for  the  Three  Flares 


Onset  Time 

Lat.  Long 

(°,o) 

V 

(km/sec) 

X 

(h) 

o 

(0) 

Case  1 

Fl:21  Sept,  04:16,  1978 

N23E40 

200 

5.0 

40 

C^ciSG  2 

F2:23  Sept,  09:44, 1978 

N35  W50 

2800 

2.0 

40 

F3:24  Sept,  17:20,  1978 

N25  W07 

m 

2.0 

50 

F6:27  Sept,  08:18,  1978 

N10W16 

550 

2.0 

50 
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Case  1 

Figure  7a  shows  both  the  observed  and  the  simulated  variations  of  ihe  solar 
wind  speed.  The  results  are  obviously  not  totally  satisfactory.  However,  this  is  the 
best  fit  obtained  by  our  method.  A  change  of  any  one  of  the  parameters  affects 
overall  results  in  very  complicated  ways,  because  the  observed  results  are  a 
superposition  of  individual  events.  For  example,  if  the  simulated  peak  speed  of 
the  solar  wind  on  29  September  is  reduced,  in  order  to  make  a  better  fit  with  the 
observation,  its  arrival  time  will  shift  from  the  observed  one,  as  well  as  the  speed 
and  arrival  time  of  all  the  other  shocks.  Thus,  the  identification  of  the  responsible 
flare  (F6)  becomes  questionable.  However,  it  was  found  that  the  other  flares  give  a 
worse  fit  than  F6.  A  better  fit  would  also  be  made  by  adjusting  the  speed  of  the 
high  speed  stream.  However,  such  a  change  affects  the  overall  variations  of  the 
solar  wind  speed  during  the  entire  period,  resulting  in  often  different 
combinations  of  the  cause  (flare )-effect  (shock)  relationship  which  is  worse  than 
what  we  have  obtained  so  far.  At  this  stage  of  development  in  this  study,  a  better 
fit  between  the  observation  and  the  simulation  is  difficult  without  additional 
observations  for  the  identification. 

The  three-dimensional  structure  of  the  expanding  shock  wave  caused  by  Fl 
is  shown  in  Figure  7b.  Note  that  the  second  shock  wave  generated  by  F2  begins  to 
appear  on  25  September;  both  the  longitude  and  latitude  lines  of  interval  5°  are 
used  to  provide  the  geometry  of  the  shock  wave.  The  shock  associated  with  the  IPS 
event  reached  the  Earth  on  25  September.  The  corresponding  IPS  sky  maps  are 
also  constructed  on  the  basis  of  Fig.  7b  and  are  shown  in  Fig.  7c.  The  crosses 
indicated  cross  points  of  the  latitude  and  longitude  lines  with  5”  intervals.  From 
Fig.  7c,  one  can  see  that  the  shock  wave  covered  almost  one  half  of  the  sky 
centered  around  the  Sun  (although  the  5^  interval  mesh  is  too  coarse  to  show  the 
shock  in  the  vicinity  of  the  center  of  the  sky  map  when  it  is  too  close  to  the  Earth). 
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Indeed,  the  IPS  observation  showed  that  this  was  indeed  the  case  (Tappin  et  al., 
1983;  a  part  of  Fig.  2  in  Hewish  et  al.  (1985J  is  reproduced  here  as  Fig.  Id.  In 
particular,  note  that  the  highest  IPS  occurred  in  the  30®-60®E  sector,  centered 
around  latitude  SO^^N.  Our  results  (Figure  7c)  are  in  good  agreement  with  the 
observed  IPS  distribution  (Figure  7d). 

Case  2 

It  is  natural  to  try  to  associate  interplanetary  events  v/ith  most  intense 
activities  on  the  Sun.  during  the  period  of  our  concern,  there  occurred  a  3B  flare 
(F2)  (N35,  W50)  at  09:44  U.T.  on  23  September;  see  Table  2.  It  has  been  generally 
believed  that  the  second  interplanetary  shock  event  (S2)  at  20:10  U.T.  on  26 
September  was  caused  by  this  3B  flare.  However,  this  association  requires  an 
average  speed  of  the  shock  wave  of  about  920  km/sec  (Cane  et  al.  1982,  1986;  Cane 
and  Stone,  1984).  One  difficulty  we  face  here  is  that  the  solar  wind  speed  at  the 
time  of  the  arr  ival  of  S2  was  only  about  400  km/sec  and  the  estimated  initial  speed 
was  more  ihan  2500  krn/scc  fCanc  et  ah,  1982).  Thus,  it  was  difficult  to  model  F2 
as  the  flare  responsible  for  S2. 

In  spite  of  such  difficulties,  we  have  attempted  to  model  the  possible 
association  between  F2  and  Si  by  choosing  the  parameters  listed  in  Table  3.  The 
corresponding  results  show  that  the  computed  IPS  disagrees  with  the  observed 
one.  In  this  case,  the  computed  IPS  was  mostly  confined  to  the  western  sky,  while 
the  observed  one  was  in  the  eastern  sky  (as  discussed  in  Case  1). 

2.2  STEP  2 

The  STEP  2  is  simply  to  compute  the  maximum  power  Pm  as  a  function  of 


V  and  B  which  are  computed  in  STEP  1.  The  power  P  thus  determined  is  a 
function  of  time. 
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FIGURE  7c 
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2.3  .STEP-ii 


(a)  Geomagnetic  Indices 

There  is  no  theory  to  relate  directly  the  power  P  of  the  solar  wind- 
magnetosphere  dvTiamo  to  the  geomagnetic  indices  AE,  Dst,  Kp,  etc.  However, 
there  are  empirical  relationships  between  P  and  AFl,  P  and  Dst,  etc. 

AE(nT)  =  -300  (log  P)^  +  11800  log  P-113200 
!  Dst  (nT)  I  =  60  (log  P-18)^  -i-  25 

Thus,  both  AE  and  Dst  can  he  computed  as  a  function  of  time  w'hen  P  is 
given  as  a  function  of  time.  Since  the  present  prediction  scheme  can  provide  Pm, 
instead  of  P,  we  can  estimate  only  the  maximum  possible  AE  and  Dst  by  replacing 
Pby  Pm. 

(b)  Cross  Polar  Cap  Potential 

ReiiT  et  al.  ( 1981)  found  that  the  polar  cap  potential  <l>pc  is  related  to  the 
Power  P  by  ■  -  (0  93P  -  319)^^^.  Thus,  it  is  possible  to  predict  as  a  function  of 
time  if  P  can  be  pr  edicted.  Since  we  can  predict  Pm,  we  shall  be  able  to  predict 
maximum  value  .  of  ^-pr. 

(c)  Equivalent  Circuit 

One  way  to  replace  the  empirical  relationship  between  P  and  AE  is  to 
consider  an  equivalent  current  circuit.  For  this  purpnr.f,  Liu  et  al.  (1988) 
developed  an  equivalent  circuit  (a  computer  code)  for  the  magnetosphere  (Figures 
8a  and  8b).  Once  P(t)  can  hi,'  inferred  from  Steps  I  and  2  (t'hgure  1),  it  is  possible  to 
determine  the  input  voltage  for  Iho  circuit.  Tlie  equivalent  circuit  can  then 

determine  all  the  magnetospheric  quantities.  Figure  8b  shows  how  the 
magnetosphere  responds  to  a  step-function  like  increase  of  The  response  as 

a  function  of  time  may  he  ;  ef'n  ch  arly  in  the  cr-.i-ss-uc-hir  cap  potentials^.,,,  and  the 

total  fc..-  rgy  deposition  hi  the  ionosphere  (denoted  by  WB). 


Polar  Cap 

Rpc 


(ujqo)  -Ly 


(A>)OOL) 


(A>)001.) 


Another  important  progress  in  this  is  that  we  have  determined  the 

relationship  between  the  cross-polar  cap  potential  drop  and  the  AE  index  on 

the  basis  of  both  satellite-based  and  ground-based  data  (Akasofu  et  al.  1990  j.  The 
results  are  shown  in  Figure  9.  Since  (hp^^.  (Kv)  -  36 0.S9  AE  or  AE  -  11.2 
(<bp(.  -36),  if  the  relationship  between  the  solar  wind  quantities  (V,  B,  0)  and  the 
cross-polar  cap  potential  can  be  established,  it  is  possible  to  predict  the  AE  index 
as  a  function  of  time. 

(e)  Polar  Ionosphere 

Once  the  cross-polar  potential  d>  (or  dJpc)  can  be  determined  as  a  function  of 
time,  it  is  now  possible  to  compute  the  electron  density  distribution  as  a  function  of 
time  over  the  entire  polar  region.  Figure  10  shows  the  computational  scheme  and 
Figure  11  shows  an  example  of  the  results.  One  can  see  clearly  that  the  ionization 
produced  by  the  solar  radiation  in  the  dayside  hemi.sphere  is  driven  into  the  polar 
cap  from  the  cusp  region.  In  the  polar  region,  however,  there  is  also  the 
ionization  caused  by  the  precipitation  of  auroral  particles.  This  mode!  reproduces 
fairly  well  the  ionospheric  trough.  Note  in  particular  that  the  trough  extends 
from  the  night  sector  to  the  day  sector,  as  discovered  by  Whalen  (1989). 

(f)  IMF  Polar  Angle  6 

An  accurate  determination  of  the  power  P  requires  not  only  V  and  B,  but 
also  the  polar  angle  8.  Unfortunaixdy,  inspite  of  our  elTorts  during  the  la.^t  several 
years,  we  are  far  from  predicting  0  as  a  function  of  time.  So  far,  our  prediction 
code  can  determine  V  and  B,  but  not  0.  Thus,  w'e  can  .‘^imply  estimate  the 
maximum  possible  pov’cr 
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In  order  to  understand  the  magnetic  field  configuration  in  the  driver  gas, 
we  have  started  to  model  a  magnetic  cloud.  Figure  12a  shows  a  study  of  the 
propagation  of  a  "magnetic  cloud"  along  a  flux  of  the  interplanetary  magnetic 
field  (IMF)  lines  (Wang  et  al.  1988;  Wei  et  al.  1990).  When  the  earth  encounters 
such  a  structure,  the  IMF  vector  at  the  earth  is  expected  to  vary  in  a  complicated 
way.  Our  study  can,  however,  predict  the  expected  type  of  variations  in  a 
magnetic  cloud  or  a  driver  gas.  It  has  been  suggested  that  such  a  "magnetic 
cloud"  is  generated  above  a  flare  and  propagate  out  in  interplanetary  space.  So 
far,  we  have  identified  several  causes  for  variations  of  0  (Figure  12b).  However, 
we  are  unable,  at  the  present  time,  to  identify  which  causes  dominate  variations 
of  0  for  a  given  flare. 

3.  PREDICTION  OF  THE  27-DAY  RECURRENT  STORMS 

3.1  Introduction 

The  past  prediction  method  of  recurrent  geomagnetic/ionospheric  storms  is 
based  on  a  study  of  the  so-called  "Bartels  Kp  musical  diagrams.”  Obviously,  such 
a  method  is  not  based  on  physics.  Like  flare-caused  storms,  the  modern  scheme 
must  be  based  on  the  prediction  of  the  power  of  the  generator.  Thus,  the 
prediction  of  the  27-day  recurrent  geomagnetic  storms  is  reduced  to  predict  the 
power  P  for  a  27-day  period  and  the  geomagnetic  indices  AE,  Dst,  Kp,  etc.  for  the 
same  period.  At  the  present  time,  it  is  not  known  how  Alfven  waves  and 
turbulence  are  generated  in  a  high  speed  stream.  Thus,  it  is  not  possible  to 
predict  the  polar  angle  0.  However,  it  is  possible  to  predict  the  of  P  =  20  VB^  sin*^ 
(0/2),  namely  I P I  =  20  VB^,  since  0  %  sin'^  (0/2)  5  1.0.  Then,  from  the  envelope  of  P, 
it  is  possible  to  infer  the  envelope  of  the  AE  index,  a  rough  trend  of  the  maximum 
AE  index  as  a  function  of  time. 
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Figure  13  shows  variations  of  the  solar  wind  quantities  during  the  period 
between  June  25  and  July  21,  1974.  Compare  the  observed  VB^  and  the  AE  index. 
We  can  see  that  the  general  trend  of  AE  index  variations  during  the  27-day  period 
follows  reasonably  well  that  of  VB^.  Note  that  fine  time  variations  of  the  AE  index 
represent  individual  substorm.  Thus,  we  can  predict  the  upper  limit  of  the  AE 
index  and  of  its  27  day  variations  if  we  succeed  in  predicting  V  and  B  for  a  27-day 
period. 

3.2  Neutral  Line  on  the  Source  Surface 

In  the  above,  we  have  demonstrated  that  the  forecasting  of  the  recurrent 
geomagnetic  activity  is  reduced  to  predicting  the  solar  wind  speed  V  and  the 
magnitude  B  of  the  IMF  for  a  27-day  period.  In  turn,  the  prediction  of  V  and  B  is 
reduced  to  inferring  geometry  of  the  neutral  line  on  the  source  surface  which  is 
an  imaginary  spherical  surface  of  radius  2.5  solar  radii.  The  natural  line  varies 
in  a  very  complicated  way  during  a  sunspot  cycle  (Figure  14a).  For  this  purpose, 
we  must  find  a  simple  way  to  reproduce  the  neutral  line.  Fortunately,  we  found 
that  the  neutral  line  on  the  source  surface  can  be  reproduced  fairly  accurately  by 
a  dipole  at  the  center  of  the  sun  and  a  few  dipoles  on  the  photosphere  (Saito,  Oki, 
Olmsted  and  Akasofu,  1989),  This  finding  has  provided  us  with  an  opportunity  to 
predict  the  geometry  of  the  neutral  line  on  the  source  surface  by  extrapolating 
time  variations  of  the  characteristics  of  the  dipoles. 

3.3  A  New  Method  of  Representing  the  Neutral  Line 

In  order  to  reproduce  the  observed  neutral  line  on  the  source  surface,  we 
assume  the  axial  dipole  at  the  center  of  the  sun  and  a  few  dipoles  near  the 
equatorial  plane  of  the  photosphere  (Figure  14b).  In  the  upper  left  diagram,  the 
observed  neutral  line  tor  Oarrington  rotation  1666  is  shown,  while  the  upper  right 
diagram  shows  its  reproduction  using  our  modeling  method  (Saito,  Oki,  Olmsted 
and  Akasofu,  1989).  In  this  particular  example,  the  two  dijxiles  are  located  in  low 
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latitudes  on  the  photosphere.  One  can  see  that  our  method  reproduces  fairly  well 
the  observed  neutral  line  in  the  upper  left  diagram.  The  lower  diagrams  are 
spherical  representations  of  the  upper  ones.  The  cross  in  the  upper  right 
diagram  shows  the  view  longitude  of  the  spherical  presentation. 


Figure  15  shows  how  accurately  we  can  reproduce  the  observed  neutral  line 
by  a  combined  axial  field  (in  the  righthand  side)  and  four  dipoles  located  in  low 
latitude,  No.  1,  2,  3,  and  4.  During  four  solar  rotations,  from  Carrington  rotation 
1661-1664,  the  four  dipoles  varied  rather  fairly  smoothly.  Figure  16  shows  time 
variation  of  the  dipole  No.  3  in  terms  of  latitude,  longitude,  the  azimuth  and  the 
magnitude.  Thus,  time  variations  of  the  equatorial  dipoles  can  be  extrapolated, 
allowing  us  to  predict  the  geometry  of  the  neutral  line,  in  turn  the  solar  wind 
speed  V  and  the  IMF  magnitude  B  and  finally  the  envelope  of  the  maximum 
variations  of  VB^  and  of  the  AE  index. 


NFUTRAL  LINE  ON  SOURCE  SURFACE 
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FIGURE  15 


The  modern  prediction  scheme  of  geomagnetic  storms  must  be  based  on 
numerical  methods,  not  on  statistical  methods.  We  believe  that  we  have 
formulated  the  entire  scheme  which  consists  of  these  computational  steps.  Our 
scheme  is  applicable  for  a  complicated  situation  in  which  many  flares  take  place 
in  a  relatively  short  period,  say  one  week.  Although  much  refinement  is  no  doubt 
needed  in  the  future,  we  believe  that  the  general  scheme  is  well  established  by  the 
present  project.  We  have  also  introduced  successfully  interplanetary’  scintillation 
obseiwations  into  the  scheme  to  increase  the  accuracy  of  the  prediction.  Most  of 
the  refinements  require  future  advance  in  the  field  of  solar  physics,  solar  wind 
physics  and  magnetospheric  physics. 

One  of  the  critical  progresses  needed  in  the  prediction  is  a  better 
understanding  of  the  nature  of  variations  of  the  polar  angle  0  of  the  interplanetary 
magnetic  field.  At  the  present  time,  it  is  not  known  how  the  variations  are 
caused,  although  some  researchers  consider  them  in  terms  of  magnetic  clouds' 
and  'driver  gas.' 

The  present  project  has  also  formulated  the  prediction  method  of  the  27-day 
recurrent  geomagnetic  storms.  We  have  shown  that  the  envelope  of  the  27-day  AE 
index  variations  can  be  predicted  by  predicting  variations  of  the  neutral  line  on 
the  source  surface  of  the  sun.  The  method  is  simple  enough  to  put  into  practice. 
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